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Abstract. We present an efficient expression for the analytic continuation to 
arbitrary complex frequencies of the complex optical and AC conductivity of a 
homogeneous superconductor with arbitrary mean free path. Knowledge of this 
quantity is fundamental in the calculation of thermodynamic potentials and dispersion 
energies involving type-I superconducting bodies. When considered for imaginary 
frequencies, our formula evaluates faster than previous schemes involving Kramers- 
Kronig transforms. A number of applications illustrates its efficiency: a simplified 
low-frequency expansion of the conductivity, the electromagnetic bulk self-energy due 
to longitudinal plasma oscillations, and the Casimir free energy of a superconducting 
cavity. 



PACS numbers: Cavity quantum electrodynamics; micromasers 42.50.Pq - Electrical 
conductivity of superconductors 74.25.Fy - BCS theory 74.20.Fg - Exchange, 
correlation, dielectric and magnetic response functions, plasmons 71. 45. Cm 
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1. Introduction 

The linear response of a liomogeneous system to an external field can be conveniently 
described by a susceptibility x(ct) ^) in Fourier space. Calculations at real frequencies 
can be directly compared to experimental spectra. Imaginary frequencies, u = i^, 
are theoretically useful to deal with equilibrium thermodynamical properties [1]. In 
particular, quantum field theories, after a Wick rotation in imaginary time, can be 
studied at nonzero temperature by expressing partition functions as Matsubara sums 
over discrete frequencies Un = 2'KniT (bosons: = 0, 1, . . .; fermions: n = |, |, . . .) [2]. 
A further advantage of this approach comes from the fact that causal response functions 
are well-behaved on the (positive) imaginary axis, a feature that facilitates both 
analytical and numerical calculations [T]. 

A strong motivation for the present paper is the interest in electromagnetic 
dispersion forces for superconducting materials. In the context of the electromagnetic 
Casimir effect, the role of absorption in a real mirror is still poorly understood. 
The superconducting phase transition provides a clean way to change absorption 
over a narrow temperature range, as has been suggested in Ref. [3]. The standard 
theory of Bardeen-Cooper-Schrieffer (BCS) ^ is well established for conventional 
superconductors and has provided the basis for calculating the optical conductivity 
a{u). In their seminal paper [5], Mattis and Bardeen included the impact of disorder 
in the 'dirty limit' where the scattering mean free path i = vpT is small. This approach 
has been generalized to arbitrary purity H E] • 

The conductivity, at imaginary frequencies, is the basic response function that 
yields electromagnetic dispersion forces between superconducting bodies at nonzero 
temperature [H [HI [ID]- In previous work, the common strategy was to compute cr(i^) 
numerically, starting from the conductivity available at real frequencies and performing 
a Kramers-Kronig transform [TT]. This technique is quite inefficient because in BCS 
theory, o"(cj) is itself an integral over the quasiparticle spectrum. This fact also 
complicates analytical studies. We present in this paper an alternative formula for 
the BCS conductivity that does not need a Kramers-Kronig transform, is given by a 
rapidly converging integral, and can be evaluated in the upper half of the complex 
frequency plane. Our derivation is based on a consistent regularization scheme where 
the quasiparticle response far away from the Fermi edge (that coincides with a normal 
conductor) is subtracted. 

The paper is organized as follows: in Sec. [2] we review the BCS conductivity at 
real frequencies. In Sec. [3] we present our formula for the BCS conductivity a{z) at 
arbitrary complex frequencies z, and we show that for z = i^, it coincides with the 
Kramers-Kronig transform, evaluated numerically from o'{u) given in Ref. [7]. In Sec. 
Hlwe discuss several applications. First, we derive the low-frequency asymptotics and 
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recover in a simple way the findings of Ref. [S]. We then compare with the two- fluid 
model of a superconductor [1] and show its limitations. As further applications, we 
study the zero temperature electromagnetic self-energy of a bulk superconductor, as a 
function of the scattering mean free path (Sec. I4.4|) . and present a numerical calculation 
of the Casimir energy between two superconducting plates near the critical temperature 
(Sec. USD. 

Our conclusions are presented in Section [51 Two Appendices close the paper: in 



Appendix A we present an alternative computation of the normal-metal conductivity. 



more in the spirit of Ref. [5]. Finally, in App. Appendix B we verify analytically that 
our formula, when considered for real frequencies, reproduces the results of Ref. [7]. 



2. Mattis-Bardeen theory 

In 1958 J. Bardeen and D.C Mattis developed a calculation of the anomalous skin effect 
in superconductors [S]. The conductivity of a bulk of superconducting material can be 
obtained from the linear current response to an applied vector potential. Expanding the 
current over the quasiparticle spectrum given by BCS theory, one finds a Chambers-like 
expression for the current density [3 H] 

- Wk R) (1) 

I{uj,R)= J de de' i^L{u, e, e') - ^^^'^ | cos[R{e' - €)/vf] (2) 

where ctq is the DC conductivity in the normal state. Disorder is taken into account 
via the scattering mean free path £, vp is the Fermi velocity, /(e) is the Fermi-Dirac 
function, and e is measured from the Fermi edge. 

The spectral function L{u,e,e') involves the quasiparticle energies {h = 1) 



E = Ve2 + A2 [A = A(T): BCS gap] (3) 
and is given by 



^E-E'-uj-iO+ E-E' + co + iO+J 
, . . J. ( l-f{E')-f{E) l-f{E')-f{E) \ 
^ 2^1^' ^ \e + E' - oj - ^ E + E' + u + ) 



with the coherence factors 

p(e,e') 1 EE'±{ee' + A'^ 



(4) 



(5) 



g(e,e') J 2EE' 

The infinitesimal parameter 0"^, set to zero at the end of the calculation, ensures an 
adiabatic switching-on of the external field with time- dependence e~^'^* and a causal 
response. The term (/ — /')/(e' ~ e) that appears in Eq.([2]) gives an imaginary 
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contribution to the conductivity (for real uj) and represents the diamagnetic (London) 
current O [12] . Note that the same expression is found by the Green function approach 
of Behtz et al. [6] when the spectral strength of the material in its normally conducting 
state is considered as diffusive within a charge- conserving approximation. 

In this paper, we focus on the local limit of the current response: it is obtained 
from the previous expression assuming that the vector potential A(r + R) is slowly 
varying in space. Taking it out of the integral and integrating over R, one obtains the 
conductivity in the form [8] 



where the scattering rate 7 = 1/r = vp/£ can be interpreted as the effective width of 
quasiparticle energies at the Fermi edge [last factor in Eq. (|6])]. This expression has been 
further simplified by Zimmermann [7] and by Berlinsky et al. |8], but their final result 
was given in a form, which involves the frequency u in the integration boundaries and in 
the argument of Heaviside functions, so that a direct continuation to complex frequencies 
is not straightforward. In the next Section, we show that the analytic continuation of 
the BCS conductivity ([6]) to arbitrary complex frequencies z in the upper half-plane C"*" 
can be written in the form 



where one recognizes in the first term the analytical continuation of the Drude result 
for the normal metal (TDr(w) = o"o/(l — iut). The BCS correction 5a-Bcs{z) is given by 
Eqs. ffTOj) . fl2Tl) below that can be easily evaluated numerically. 

3. Analytic continuation to complex frequencies 

Let us now go back to the expression given in Eq. We note that the spectral 

function (j4]) is evidently analytic for u = z & C"*", and the analytic continuation of the 
BCS conductivity ([6]) is immediate at this stage. The double integral over e and e' is, 
however, tricky to perform. 

Mattis and Bardeen rewrite the integrand using symmetries under the exchange 
of e -f-)- e' and introduce a regulating function to ensure the convergence at large e, e' . 
This makes obsolete the inclusion of the diamagnetic term under the integral in Eq. ([6]) . 
They then perform the integration along a contour in the complex e'-plane, leaving only 
the integration over e to be done numerically. Further calculations are performed in the 
'dirty' limit where the mean free path i is much smaller than the coherence length w^/A, 
i.e., 7 ^ A. Zimmermann [7] and Berlinsky et al. [S] consider a general scattering rate, 
but focus on the case of real frequencies. 





(7) 




(8) 
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In the following, we show that the difficulty with the convergence can be overcome 
by a proper handling of the integration scheme: we exploit the exchange symmetry 
in the ee'-plane and secure uniform convergence of the double integral by a suitable 
subtraction. We take into account explicitly the diamagnetic current and include a 
general scattering rate as in Refs. [71 18] . This procedure makes the Drude term in Eq. ([8]) 
emerge in a natural way for complex frequencies and leads to a numerically convenient 
form for the remaining BCS correction (whose integrand is strongly localized). 



3.1. Normal metal 

Let us consider first the limiting case of a normal conductor, that is the expression 
obtained from Eq.(l6]) in the limit A — 0. We allow for complex cj, denoted by z. The 
coherence factors in Eq.(|5]) reduce to |[1 ± sgn(ee')]. We observe that the first and 
second terms of Eq.(j4]) are related by the mapping E' -H- —E'. It is easy to check that 
in all four quadrants of the ee'-plane, we can replace the spectral function by 

A ^ : L{z, e, e') ^ [/(^) " /(^')] (9) 

The conductivity of the normal metal is now obtained as follows. Introducing the 
variables s = |(e + e') and 5 = e' — e, the difference of Fermi functions can be written 
as [ks = 1] 

m - m = . ,.;"'':'f ^' „. (10) 

cosh(d/2T) + cosh(s/T) 
which is integrable at \s\ — t- oo with the result 

H-oo 

/d.[/(e)-/(e')] = 5 (11) 

— oo 

The integration over 5 is then easily performed: 

.W = !5!l7dW-^-l|B,W = ^ (12) 

iz J \o'^ — z'^ I 7 — 12; 

—00 ^ ^ 

where the subtraction in curly brackets corresponds to the limit 2 — (the diamagnetic 
current, compare Eqs. ([6l [9])). We thus find the analytic continuation of the Drude 
conductivity. Appendix A| presents an alternative derivation of this result, that proceeds 



in close analogy to Ref. |5]. In this approach, the integrand is multiplied by a factor 
exp[i(e + e')/(2A)] in order to enforce absolute convergence in the direction |e + e'| — 00 
and the limit A — 00 is taken in the end. 



3.2. Superconductor 

We now return to the full BCS theory. It is a simple matter to verify that for any 
complex frequency z G C"*", the integrand involving the spectral function L(2;,e, e') can 
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be split into two terms 



L{z, e, e')D,{t - e') = F{z, e, e') + F{z, e', e) , 



(13) 



where 



F{z,e,e') 



D,{e-e') 



tanh{E/2T) j A+ + €€' + ee' 



(14) 

(15) 
(16) 



A±{z,E) 
Ql{z,E) 



E{E ± ^) + 
{E±zf-^^. 



Eq.( !T3|) shows manifestly the symmetry of the integrand under the exchange e O e'. 
Mattis and Bardeen proceed by integrating separately over the two terms F{z, e, e') and 
F(z, e',e), but to do so, one must ensure absolute convergence of the double integral. 
Indeed, for large |e|, |e'|, the tanh functions tend to ±1, and from Eqs. (11 411TB]) . we 
see that F(2,e,e') becomes proportional to (e — e')~^. Therefore for large |e|, |e'|, the 
integrand F[z,e,e') becomes independent of the variable s = (e + e')/2, and then its 
double integral is not absolutely convergent as |s| — )■ oo. 

This divergence can be conveniently cured by the following subtraction procedure. 
We evaluate the asymptotic limit of F{z,e,e') as |e + e'| ^ A, keeping the ratio to 
the other parameters z, T, 7 fixed, and subtract it from F. This yields a regularized 
integrand 



We have checked that this difference is of order 0[(|e| + |e'|)~^] in the direction e = e' 
and 0[(|e| + |e'|)~^] in all other directions of the ee'-plane. In addition, there are no 
poles for real-valued e, e' so that the integral of F{z, e, e') is uniformly convergent. 

Comparison of Eqs.(j9l [17]) shows that this procedure is exactly equivalent to 
subtracting the spectral function of the normal conductor from the BCS expression 
L{z, e, e'). We hence find a form of the conductivity as anticipated in Eq.(]8]), where the 
BCS correction is 



In the second line, we have exploited the absolute convergence of the integral of 
F[z, e, e'), to interchange the order of the e , e' integrations for the second term. 

The BCS correction Sascsi^) [see Eq.(jH])] is now evaluated as has been done by 
Mattis and Bardeen for real frequencies. We remark however that since the subtracted 
integrand in toto converges uniformly, no cutoff is required. The e' integration of the 




(17) 




(18) 
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quantity F{z, e, e') is easily performed with contour techniques. We first note that thanks 
to the subtraction, the integrand is at most of order |e'|~^ as |e'| — )■ oo at fixed e so that 
we can close the contour by a semi-circle at infinity in the upper half-plane. We then 
observe that the second term between the curly brackets in Eq. (fT7|l is an odd function 
of e' — e and vanishes upon integrating over e' . The first term has simple poles occurring 
at e' = e + 17, and at e' = Q±. The square root Q± must be taken with a cut along the 
real positive axis [13], in such a way that Im(Q±) > 0. Evaluation of the corresponding 
residues gives the result: 

5(Tbcs{z) = i^ ^ tanh (^^^ ^ ^^(z, e) , (19) 



r,=± 



where 



We note that both terms in Eq.fl20l) have a pole at e = Q± — 17, but the residues 
compensate each other exactly. This creates numerical problems, that can be avoided 
by rearranging into 

. e^Q±iz,E) + iQ^iz,E)+i^)A^iz,E) 
' ^ QAz,E)[e^-{Q^{z,E) + i^y] 

where we have restored the arguments of Q± and A± for clarity and E = E{e) as in 
Eq.03]). Note that the integrand in Eq. f|T9|) is actually even in e. 

Eqs. (fT9| |2T]) . together with Eq.([8]), provide the desired analytic continuation of the 
BCS conductivity to complex frequencies z G C"*". In Appendix B we verify that our 
formula, when considered for real frequencies u, correctly reproduces the expression for 
cr(ci;) quoted in Refs. [3 [8]. 

Consider now the special case of a purely imaginary frequency z = i^, that one 
needs for thermal equilibrium quantities calculated in the Matsubara imaginary-time 
formahsm. By using the relations Q-{i^, e) = — (5+(i^, — e), and A_{i^, e) = A*^{i^, — e), 
it is easy to verify that the functions G±{iC,, e) satisfy the relation G^{i^, e) = G^(i,^, — e). 
Then, from Eqs.f ll9p and (120 p we obtain: 



SaBCsm = / § tanh(E/2T) Re[G+(i^, e)] (22) 

—00 

In the normal conductor limit (A — 0), G+{i^,e) becomes purely imaginary so that 
Eq. (l22!) vanishes, leaving only the Drude term in Eq.([8]). 

From this expression we can obtain the 'extremely anomalous' or 'dirty' limit 7 ^ A 
considered by Mattis and Bardeen, by expanding in powers of I/7, and retaining only 
the lowest order: 

7<3+(i^,e) 
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Figure 1. (Color online) BCS conductivity evaluated along the imaginary positive 
frequency axis. Lines: Direct analytical continuation derived in Eqs.Q and (|2H 
[221) ill the clean (7 « 1.4 A, black line) and dirty (7 « 1400 A, gray line) case and 
the extremely anomalous limit (7 ^ A, red line) using Eg. ([25)1 . Dots: Kramers- 
Kronig continuation from the real frequency conductivities by Zimmermann et al. [7] , 
Bcrlinsky et al. [S], following the scheme summarized in Eq. ([M)) for the clean and 
dirty case. Diamonds: Kramers-Kronig continuation in the extremely anomalous case 
starting from the complex conductivity (t(w) of Ref. 



In this limit, the BCS correction (Jctbcs [Eq. fl22l) ] no longer depends on 7. Since the 
approximation is valid for small frequencies ^ ^ 7, the Drude contribution in 
becomes (TDr(iO ~ '^0 in this limit. 

3.3. Numerical efficiency 

Kramers-Kronig relations are commonly used to calculate conductivities or dielectric 
functions at imaginary frequencies starting from the known representations at real 
frequencies (see, e.g., Ref. [H]). This is particularly useful if one wants to extrapolate 
experimental data to the imaginary axis. The formulation we use here. 



involves the imaginary part of the conductivity (an odd function in u). We thus avoid the 
6{u) singularity [5l[8] that occurs in Rea{u). Eq. (!2^ is numerically quite inconvenient 
in the BCS theory since it requires two integrations (over e and u). The calculations are 
more easy for the direct analytical continuation, Eqs.(l8| [2T | [221) . where a single integral 
must be done that in addition converges rapidly. The agreement between the two 
representations is excellent as one can see from Fig. [1] The solid lines are obtained with 
our analytical continuation while the dots result from the Kramers-Kronig scheme. Our 
results also reproduce the dirty limit based on Eq.(l23l) (7 ^ A) considered by Mattis 
and Bardeen [5] for real frequencies. All the numerical work has been performed using 








(24) 
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the approximate form of the BCS gap function A(T) given in Refs. [TH [T5] . 

4. Applications 

4.1. Supercurrent response 

As a first example, we show that the BCS response function can be analyzed in a quite 
simple way at imaginary frequencies. We recall the discussion of Berlinsky et al. [8] 
who attribute the weight of the peak oc 5{uj) in Recr(w) to the Cooper pair condensate 
and the response at nonzero frequencies to thermally excited quasiparticles [u <^ 2A) 
and to 'direct' excitations across the gap {u > 2A). All three contributions together 
fulfill the sum rule for the oscillator strength in the conductivity. The first and second 
contributions are associated to the superfiuid and the normal fluid of the two- fluid model 
due to Gorter and Casimir [161 S]) whose weights vary with T in a way similar to the 
Gorter-Casimir rule. The third contribution is found to be approximately constant over 
a large range of T. 

We have tried to evaluate the imaginary frequency representation of these three 
contributions, taken separately. We recognized, however, that the 'thermal' and 'direct' 
parts of the quasiparticle response [i.e., first and second term in Eq.(j4])] show a branch 
cut in the e'-plane so that only their sum can be combined in the same way as in 
Eqs. ffTBHTB]) . In what follows, we hence consider a slightly different 'two-fluid split' of 
the conductivity that is closer to the Gorter and Casimir approach: a supercurrent 
part that does not show any relaxation (or scattering) and a resistive current whose 
frequency response is Drude-like. In addition, we reproduce the T-dependent weight 
of a logarithmically divergent contribution to the conductivity at low frequencies given 
in Ref. [8]. This behaviour is usually attributed to the divergent mode density of the 
BCS approach at the gap that is not smoothed out in the Chambers-like way disorder 
is introduced by Mattis and Bardeen [5]. 

Along the imaginary frequency axis, the superfiuid response corresponds to the 
weight of the pole at ,^ = 0. In terms of an effective superfiuid plasma frequency. 



where the logarithmic term is discussed below. 

Comparing to Eq.(IH]), we get Us by considering (5crBcs(iO the limit ^ — 0. Our 
starting point is then Eq.( l22l) . For { — )■ 0, the integrand G+(iC, e) [Eq. (l2T]l ] as a function 
of e exhibits two relevant features. First, a 'narrow peak' originating from the term 
l/Q+(i^,E(e)): its width in e scales like (^A)^/^^ Setting e = ^(^A)^/^ we get 



+ i?(T) log (A/0 + C(T) 



(25) 




2f7oA 



/ 



dx tanh 



Re 



2i 



5^0 o"oA 



tanh(A/2T) (26) 



2T 



Vx^ + 2i 



TT 







Dispersion energies of BCS superconductors 



10 



where convergence at the upper hmit is secured by taking the real part. In the last step, 



we have taken the lowest order in the expansion of E{e) = Ay 1 + (^/A)x^ for small ^. 
In the limit T — )■ 0, we recognize here the integral term in Eq.(14) of Berlinsky et a/. [8]. 

The second feature in ^^(i^, e) is a 'broad background' whose width as a function 
of e is set by A and 7. This gives a negative contribution to u^: 

tanh(E/2r) 



CO. 



s,2 



4^A2 

4^A(0) 

£0 







VA2 + e2 (72 + 4e2) 
arcsec(2A(0)/7) 



(27) 



/4-(7/A(0))2 

The zero-temperature limit given in Eq.(15) of Berlinsky et al. [8] is identical to the sum 
of Eqs.( pBP7|) . Recall that in this paper, A always denotes the temperature-dependent 
gap. 
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Figure 2. (Color online) Superconducting order parameter for different disorder 
strengths. We plot the superfluid part of the current response function, expressed 
in terms of the superfluid plasma frequency squared, w^, and normalized to the 
weight of the high-frequency (above gap) response = (Tot/eo- This is plotted 
as a function of the temperature (top) and of the scattering rate 7 (bottom). Dashed 
line: normalized BCS gap A(r)/A(0). The superconducting fraction vanishes above 
the critical temperature ~ 0.56 A(0) and becomes almost constant for T ^ Tc- In 
this limit, ojg still depends on 7 and reaches the normal-phase plasma frequency tu'^i 
only in the pure limit 7 — > 0. 



Figj2] shows the behavior of = + UJI2 as a function of temperature and 
scattering rate. We normalize to the plasma frequency given by the total carrier 
density, cUpj = cto7/£o- This suggests an interpretation in terms of a superfluid order 
parameter [161 [H] • Indeed, the superfluid response identically vanishes for temperature 
higher than the critical temperature and goes to a constant in the limit T <^ Tc. 
Note that at low temperatures, Ug still depends on the disorder in the sample via the 
parameter 7 (see also in the following). 
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4-2. Logarithmic correction 

It is well known that at nonzero temperature, the BCS conductivity shows a 
logarithmically divergent term at low frequencies. The weight B{T) of this term [see 
Eq. (!25l) ] can also be calculated along imaginary frequencies, by improving on the small- 
^ limit discussed before. The 'narrow peak' of Eq. fl26|) allows to take into account the 
small e expansion of the temperature-dependent factor tanh£'/2T. The 0{e^) term of 
this expansion gives then the following integral in the scaled energy x = e/(^A)^/^: 



T 



Xc 


ix^ 


J dx Re 




Vx2 + 2i_ 



(2^ 



A cutoff Xc is needed for the integral in order to reproduce the convergent behavior of the 
original expression. Inspection of the full integrand shows that it starts to decrease at 
e > A so that we take Xc = (A/^)^/^. The integral can then be evaluated explicitly and 
gives to leading order at small ^, logXc = | log (A/^). We thus read off the coefficient 
B{T) in Eq.(I25]): 

(XoA 



B{T) 



2T 



■sech^(A/2T) 



(29) 



It is easy to check that in the regime u ^ T <^ A, our result (analytically continued to 
real frequencies) overlaps with the formulas (20, 21) of Berlinsky et al. [S]. Their result 
apparently allows for larger frequencies up to the order T, but is restricted to w, T ^ A. 

4-3. Partial sum rule 





Figure 3. (Color online) Effective oscillator strength in the frequency range < cj < ^, 
as measured by the ratio ^cr(i^)/(cro7)- The thick red curve is based on the full BCS 
conductivity, the thin black one on the Drude part only, ^/(7 + ^), and the thin red 
one on the BCS correction. Dash-dotted line: two-fluid-like model of Eq.(|30|) with the 
'order parameter' t]{T) matched to reproduce the small-^ limit. 

We consider a 'clean' and a 'dirty' case at a low temperature (top) and at a higher 
temperature (bottom). Note that A is always the T-dependent gap here. The shoulder 
at low frequencies in the 'clean and warm' case (right, top curve, 7 = 0.2 A, T = 0.45 A) 
is probably related to the logarithmic contribution in Eq.p5|). 
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We now make a connection between the conductivity at imaginary frequencies and 
an effective oscillator strength, integrated over in a finite range of frequencies. More 
precisely, we consider the product ^o"(i^) that is plotted in Fig. [31 A Kramers-Kronig 
relation, similar to Eq.( l24l) . but based on Re o"(co'), shows that ^o"(i^) is an non-decreasing 
function of ^ that collects the oscillator strength in the frequency range < w < The 
sum rule for Recr(a;) thus implies that ^o"(i^) — )■ 7cro for ^ — oo. According to Eq. (!25|) . 
the limiting value of Ccr(i^) for ,^ — )■ is SQuliT) =: r]{T)jao where < r]{T) < 1 
is similar to a superconducting order parameter. A normal conductor necessarily has 
r](T) = because its DC conductivity is finite. 

We compare the product ^o"(i^) in Fig. [3] to an effective two-fluid model with the 
conductivity (dash-dotted lines) 

c^2F(iO = r/(T)^ + (1 - r/(T))^ (30) 

where the damping rate for the normal fluid contribution (second term) is set to 7. The 
overall agreement is quite good, but the logarithmic correction is of course not captured 
by this approach. We also note that the 'order parameter' t]{T) introduced in this way 
does not coincide with the gap function A = A(T), see Eqs.fl26[ [27]) . In particular, it 
depends both on temperature and the dissipation rate and does not reach ri{T — > 0) = 1 
as long as 7 is nonzero (see also Fig. |2]). 

4.4- Plasmon self-energy 

As a further example, we consider the electromagnetic self-energy of a bulk 
superconductor at zero temperature. This quantity can be calculated from the oscillator 
strength of the bulk plasmon oscillation [18]. Recall that plasmons are collective, 
longitudinal oscillations of the carrier density that are up-shifted in energy through 
the Coulomb interaction [191 ISO]- We calculate here the self-energy in the local limit 
and find the longitudinal (density) response from the dielectric function 

.(c.) = l + i^ (31) 
Sou 

The bulk plasmon dispersion relation is indeed given by s{(jj) = with solutions in 
the lower half-plane C^. Along the real frequency axis, the spectral density broadens 
and shifts, and applying the logarithmic argument theorem, the self-energy is given 
by [m mi [221 123] 

Performing a partial integration and shifting the integration path to imaginary 
frequencies, we get 

r°° duj f°° d^ 



E 



Im/ — log e{uj)= — loge(iO (33) 
JO Ztt Jo ivr 
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where the log^ divergence at ^ — ?■ is integrable. For a normal conductor [e{uj) in 
Drude form], this is easily evaluated in terms of two characteristic frequencies, f2pi and 
7: 

^Pi = V'^pi-7V4-i7/2 (34a) 

where Wpi = ^J'yao/eo is again the plasma frequency. The case 7 <^ Up\ is the relevant 
one for most materials. In the limit 7 — 0, we recover the zero-point energy U!p\/2 
for an undamped plasmon [dielectric function e{uj) = 1 — 00^1/ u'^]. At nonzero 7, 
corrections are due to the renormalization of the real part of Qpi and the logarithmic 
term. The latter can be interpreted as the contribution of zero-point fluctuations of the 
bath that is responsible for the damping of the plasmon mode [25] . Alternatively, 
it represents how the modes of the bath shift in frequency, due to the coupling to 
the plasma oscillator [18]. We quote here only for completeness the overdamped limit 
(wpi < 7/2) 

ED = -^^\og^j, e± = 7/2 ± V7V4^ , (346) 

although the Mattis-Bardeen theory for an impure superconductor [S] is clearly no longer 
valid in this limit. 

For the superconductor, the self-energy E = Es follows very closely the normal 
conductor, as can be seen from Fig. HI The difference arises from collective modes within 
the gap and modifies Es by an amount of order A. The plot features the transition to 
the dirty limit [5] at 7 ~ A(0) and a saturation of Es — E^ in the overdamped limit. 

A 

A 
2 



0.1 1 10 100 1000 

J / A 

Figure 4. (Color online) Electromagnetic self-energy Es [Eg. ([33| with the BCS 
dielectric function] of the bulk plasmon mode in a superconductor vs the impurity 
parameter 7 (logarithmic scale), at zero temperature. Thick black line: self-energy 
Es, in units of the plasma frequency LUp\ « 2800 A (0) (left scale). Red line: difference 
Es — Eo to the normal conductor, Ea. (|34ap . in units of the (zero temperature) BCS 
gap A = A(0) (right scale). We take for Ed a Drude dielectric function with the same 
plasma frequency Wpi and damping rate 7. 
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D 



Re 



- -pi . - -pi 1 - -pi 

— -1 — log^ 

2 TT 7 
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4-5. Casimir free energy of a superconducting cavity 

Another example where the knowledge of optical response functions at imaginary 
frequencies is not only highly useful, but virtually indispensable, are electromagnetic 
self-energies like the Casimir interaction. Indeed, the calculation of these energies is 
practically impossible to perform along real frequencies because of a highly oscillatory 
integrand. Along the imaginary frequency axis, the integrand becomes smooth and 
rapidly convergent, similar to what happens in Eq.( l33l) . But the Casimir effect in a 
superconducting cavity is very promising also from a physical viewpoint because the 
superconductor offers the possibility to control dissipation. The influence of the latter 
on the Casimir free energy remains indeed an open question [101 EZ] : two conflicting 
viewpoints have been raised as to whether the DC conductivity of the mirrors should be 
included in the modelling or not. It has therefore been proposed to resolve this issue by 
considering the Casimir energy of a superconducting cavity, more precisely its change 
across the critical temperature [3l[TTl[28]. Earlier calculations in the framework of BCS 
theory [11] had to use a Kramers-Kronig analytic continuation of the optical response 
functions to imaginary frequencies [Eq. flM]) ]. leading to a significant computational 
overhead. The scheme developed in this work involves requires less integrations and 
enables us to perform numerical analysis in a much more flexible and precise way. 

We recall that the Casimir free energy per unit area [9] between two plates 
(separation L, temperature T) is given by a sum over the Matsubara frequencies 
= 27rnT, 

^(L, T) = — ^ 7 kdk log Dj,{tU k) (35) 
zvr JO p 

Dp{uj, k) = 1 — rp{u, k) exp (— 2kL) , 

where the primed sum weighs the zeroth term by 1/2 and the last sum is over 
polarizations p G {TE, TM}. The simplest model for the reflection amplitudes are 
the Fresnel equations 

rTE(c^, k) = — ■ , rTM(w, k) = — - — ■ (36) 

K + Km e{U)K+Km 

with the notation k = ^/k"^ — ^ and = \Jk'^ — ^i^)^ (for ^ £ C'^, Rck > and 
Im/t < 0). In the numerical calculation of the Matsubara sum from Eq. (j35|) . the sum 
has been cut off at a sufficiently high [29] value of n and the zeroth summand (^o = 0) 
has been evaluated using the results of Sec. 14.11 

Fig. |5] shows the temperature dependence of the Casimir free energy calculated 
for different theoretical descriptions using the parameters for niobium. Both the BCS 
description and the simple two-fluid model [Eq. llSU]) with ri(T) = 1 — (T/Tc)^] predict 
a sudden change of the Casimir energy across the superconducting transition. As 
estimated in Refs. [5], [HI |2B] , this effect should be observable in experiments. 
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Figure 5. Casimir free energy (per unit area) vs. temperature in a planar cavity of 
length L = 100 Ap where Ap = 2nc/uj-p\ is the plasma penetration depth. For numerical 
work we used values for niobum [3D] i^pi = 1.408 x lO^^rad/s, 7 = 2.44 x lO^^'rad/s, 
Tc = 9.25 K, and A(0) = 1.7 « 8.5 x 10"''7 (dirty limit). Energies are normalized 
to ~ 0.2Tc/L^, the absolute value of the free energy in the Drude model at T^- 
The Matsubara sum ([35]) is computed including up to 10,000 terms at the lowest 
temperatures. 



Having in mind the discussion from Sec J4.1l it is not surprising that there is an offset 
between the BCS curve and the two-fluid model (or plasma) curve at low temperatures. 
This is because impurity scattering renormalizes the supercurrent plasma frequency cOs, 
leading to a smaller effective oscillator strength in the BCS conductivity, in particular 
for low ^. 

5. Conclusions 

In this article we have analyzed the analytical continuation of the bulk BCS conductivity 
to imaginary frequencies. The final result is numerically more efficient than previously 
used schemes based on Kramers-Kronig relations. Our approach simplifies the direct 
calculation of most of the involved integrals, and displays clearly convergence issues. It 
also illustrates how the conductivity of the normal component emerges naturally from 
the BCS framework. The superconducting correction vanishes as the BCS gap A — > 
and is given in terms of a rapidly convergent integral. The result is convenient for 
both analytical and numerical computations of the optical response of superconductors. 
We have illustrated this by re-deriving in a simple way the weight of the supercurrent 
response, including its dependence on temperature T and the mean free path i = vp/'j, 
and the weight of the logarithmically divergent term in the conductivity. As another 
illustration, we have given numerical calculations of the change in the bulk plasmon 
self-energy at zero temperature, showing how impurity scattering increases the self- 
energy, although the superconducting gap lies well below the bulk plasmon resonance. 
Finally we have demonstrated how the Casimir energy for the case of a cavity made of 
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superconducting material can be calculated numerically in an efficient way and recovered 
the energy change in the superconducting transition. 

A more detailed investigation including the discussion of the distance dependence of 
the free energy and possibilities to tune Casimir energies through the dissipation rate 7 
will be presented elsewhere. Future work will also discuss the temperature-dependence 
of the self-energy. Preliminary studies have revealed features (related to the behaviour 
of o"(i^) near = 0) that lead to similar physical effects as in the thermal correction to 
the Casimir energy between metals or ideal conductors [lOl |26l |27]. It would be also 
interesting to extend the calculations beyond the local (or macroscopic) limit, taking 
into account a finite wave- vector q, as done by Popel [31]. This would also provide 
a natural scheme for calculating the electromagnetic self-energy of a superconducting 
surface [321 [33] • 
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Appendix A. Alternative calculation for the normal metal 

This approach is closer to the way in which Mattis and Bardeen performed the 
integration. We start by combining Eqs.(|6l |9]) and get for a complex frequency z G C'^: 

'"dede' [tanh(e72r) - tanh(e/2r)] 



"^^^ = 2i7 



X 



D^(e-e')e'(^+^')/2^ (A.l) 



where the limit — of the integral has still to be subtracted. We have introduced a 
convergence factor e^'y'^+'^')l'^^ which ensures absolute convergence also in the sum variable 
e + e', when analytically continued into the upper half plane. The order of integrations 
is now immaterial. Following Mattis and Bardeen, we can handle separately the two 
terms with the hyperbolic functions: integrate the term proportional to tanh(e'/2T) 
over e (closing the contour in the upper half of the complex e-plane) and over e' the 
term involving tanh(e/2T). Summing the two results, this gives 

~a{z) = ^ /detanh(e/2r)^^^ (e-/^^ - e-'^') (A.2) 

Due to the convergence factor, the integral over e does not vanish by parity. It can be 
calculated by considering a rectangular contour in the complex e-plane that encloses a 
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strip of height 27rT in the upper half-plane. In this way, we find 

cro7 27riT iri'-f 



(Tiz 



niz sinh(7rT/A) z'^ 



-7/2A _ iz/2A 



r 



(To7 



1^(7 



iz 



(A.3) 



where the limit A — )■ 00 was taken in the last step. We still have to subtract the limiting 
case 2; — )■ 0, i.e. the diamagnetic term. This gives the Drude conductivity CDrl^;) of 
Eq.(Il2D. 



Appendix B. Comparison with Zimmermann's formula 

In this Appendix, we prove that our formula for the conductivity, Eqs.(|H]), (IT^ and 
f l2U]) . when considered for real (positive) frequencies u is equivalent to the following 
expression, quoted in Ref. [7], and reproduced here for the convenience of the reader 
fwe take ^ = 1): 



criu) 



2uj 

J{uj < 2A) 
J{uj > 2A) : 



i^(J 



tanh ( —=] 
\2TJ 



dEh , 
dEh + 



dEh 



(B.l) 

(B.2) 
(B.3) 



17 



17 



, (B.4) 



tanh 



E 



2T 



A, 



P1P2 



P1-P2 



17 



A, 



P1P2 



+ tanh (A) 
= tanh (|;) 



A, 



P1P2 

A_ 



Pi 



P3P2 



P2 
1 



17 



1 + 



A, 



P3 + P2 + i7 



P1P2 

A_ 



-Pi 
1 



P, 



17 



(B.5) 



Pi 



P3P2 



^2 + 17 , 
1 



P3-P2+h 



,(B.6) 



where 



Pi 

P3 



(E + u) 



A2 



00 



- A2 



P2 
Pa 



_ A2 



iVA2- (E-w) 



(B.7) 
(B.8) 



In the relevant ranges of w, the arguments of the square roots are non- negative real 
numbers, and the square roots are then defined to be positive. This conductivity differs 
from our Eqs. ([8]), (fT9l) and (I2T1) by the presence of the frequency u in some of the 
integration boundaries. To demonstrate the equivalence with our formulae, we show as 
a first step that Eqs. (lB.ltlB.8P can be split into a Drude term and a BCS correction, as 
in Eq.®. We note that the quantities Pi, P3 and P4 above are related to our quantities 
Q±{ijj,E) as follows: 

Pi(a;,E) = g+(u; + iO+,E) , (B.9) 
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Psiuj > 2A, E) = Q^{uj + iO^,E) for A < E < u 



A 



and 



P^iuj < 2A,E) = g_(w + iO+,E) for A < ^ < w - A. 
Piiu > 2A,E) = g„(w + iO+,E) hi Lo - A < E < uj + A. 
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(B.IO) 

(B.ll) 
(B.12) 



Using relations ( ]B.9tlB.12|) above, we verify that for all non-negative real frequencies u, 
Eqs.f lB.lllB.SP can be recast in the following compact form 

icro7 



where 



aiuj 



2cu 



-H(u) 



dE 



tanh 



5a{uj) , 
E + u 



(B.13) 



2T 



B_{uj,E 
E \ 



and 



with 



-^(i?-A-u;)tanh(— j B^oo^E) 



^ f E \ 



(B.14) 
(B.15) 



a=± 



B±ico,E) 



A, 



Q±P2j Q± + P2 + h V Q±P2j Q±-P2 + h 

We can now show that the first term on the r.h.s. of Eq. flB.13l) coincides with the 
normal metal contribution. For this purpose, consider the quantity H{uj). If the integral 
of tanh(£'/2T)i?_((jj, -E) were absolutely convergent at infinity, it would be possible to 
perform the shift E + uj ^ E m the integration variable E in the first term between the 
square brackets of Eq. flB.14p . and then we would find that H{(jj) is zero. However, since 
for large values of E, the quantity taii]i{E/2T)B^{uj,E) has the asymptotic expansion 

2 



(B.16) 



tanh(E/2T) B^{u,E) 



0{E- 



(B.17) 



a; + 17 

its integral diverges at infinity. Therefore the shift of integration variable is not 
permitted, and the conclusion that H{uj) vanishes is not warranted. To evaluate H{u), 
consider the following cut-off function fh{E): 

A2 



E^ + A^ ' 



(B.18) 



Thanks to the difference appearing in the square brackets of Eq.( ]B.14l) . the integral over 
E is absolutely convergent at infinity, we have the identity: 



H{oj) = lim / dE 



tanh 



E + u 
2T 



-e{E-A-uj)tanh{ — )B4oo,E) 



E 



B^{uj,E + uj) 

jAiE) . 



(B.19) 
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Since, for any A < oo, the cut-off integrals of the two terms between the square brackets 
are individually convergent at infinity, it is now legitimate to shift the integration 
variable in the first term. After we do that, we get: 

H{uj)=\im r dEtanh(^)B_{co,E)[U{E-u)-U{E)]. (B.20) 

For large values of A the difference between the cut-off functions is appreciably different 
from zero only for large values of E, and therefore, we can replace tanh(£'/2T)i?_(a;, E) 
by its large- expansion. Eg. (IB .17]) . Inserting into Eg. (IB. 20]) . the evaluation of the 
elementary ^'-integral results into: 

Hiuj) = (B.21) 

After substituting this expression for H{uj) into Eg. (]B.13l) . we indeed find as promised 
that its first term reproduces the normal-metal conductivity (707/(7 ~ i'^)- This shows 
that Egs. (IB.ltlB.Sp for the BCS conductivity a{uj), once recast in the form of Eg. (lB.13p . 
are indeed analogous to our split form Eg. ^ of the conductivity. 

What then remains to show is that our expression for 5a-Bcs{.z), given in Egs. lIT^ . 
(l^T]) . when considered for real freguencies, coincides with the guantity 5a{(jj) in 
Egs. OB.lSp . (lB.16p . This can be worked our in a straightforward way by making the 
change of variables e — sgn(e)-E in Eg. flT^ . We omit details of this lengthy calculation 
for brevity. 
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